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Abstract 

We construct a rigorous renormalisation scheme for analytic vector 
fields on T 2 of Poincare type. We show that iterating this procedure 
there is convergence to a limit set with a "Gauss map" dynamics on 
it, related to the continued fraction expansion of the slope of the fre- 
quencies. This is valid for diophantine frequency vectors. 
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1 Introduction 



A renormalisation operator 1Z is defined for analytic vector fields on the 
torus of dimension two, associated with a frequency vector. This approach 



is based in |T(J for Hamiltonian functions, and in \12\ for flows on T d = 
M. d /Z, d , d > 2. The vector fields considered here are of Poincare type, i.e. 
there is a classification by a unique winding ratio on which 1Z acts as the 
Gauss map. The slope a of the winding ratio is mapped by the shift of its 
continued fraction expansion. All vector fields corresponding to non-zero a 
are renormalisable, and those with irrational a are infinitely renormalisable. 



We extend a result in [12] based on [1C] for diophantine winding ratios. 



The renormalisation group technique has been used for some time ap- 
plied to various problems in dynamical systems (e.g. |7], 13, 17|). Roughly, 



it analyses systems on longer time scales and smaller spacial scales, gener- 
ating a new system. By iterating this we may get convergence to a limiting 
behaviour, preferably simple, e.g. a fixed point. In such a case, we get 
self-similarity, which can be non-trivial depending on the problem. 

We show here that an orbit by 1Z of a constant vector field with dio- 
phantine winding ratio, attracts all the nearby orbits in the same homotopy 
class. This can be applied to other flows on domains with an extra vertical 
dimension. Also, it might be used for two-degrees of freedom Hamiltonian 
systems, in order to extend p0| ]. The allowed set of frequencies is smaller 
than the one obtained by KAM theory. Nevertheless, this is not a real 
disadvantage of our method since it is of equal full Lebesgue measure. 



Our procedure can be related to the one followed by MacKay [14], and 
Chandre and Moussa || for an approximate renormalisation scheme in the 
framework of the breakup of invariant tori in Hamiltonian systems. The 
renormalisation studied by those authors, besides truncating the Fourier 
modes of the Hamiltonians, does not consider perturbations depending on 
the action variable. That would have implied a change in the Hamiltonian 
vector fields, namely the angle time derivative. Thus, it does not have the 
unstable direction that appears in our systems corresponding to perturba- 
tions non-collinear to to. 

In Section [||, we include the main idea for the construction of the linear 
part of the renormalisation, along with the main properties of the continued 
fraction expansion of an irrational number. Then, in Section ||| we introduce 
the space of vector fields of our interest. We define a one-step operator 1Z in 
Section |I| and show one of its statements in Section 0. It follows Section || 
with more known results of number theory, this time on the diophantine con- 
dition. The main convergence theorem for diophantine vectors is presented 
in Section ||, along with a proof in Section |8[ 
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2 Change of Basis for a Generic Frequency and 
the Continued Fraction Expansion 

It was shown in [12j (using results from jOJ) that for a quadratic irrational 
slope of a two-dimensional vector uj, the renormalisation operator therein 
constructed is well-defined since its linear part exists. More precisely, we can 
find a hyperbolic matrix T in GL(2,Z) for which u is an eigenvector. For 
other vectors one needs to construct a different operator. In fact, for each 
iteration of the renormalisation we will consider a different transformation, 
determined by the continued fraction expansion of a. 

Consider the constant vector field u on the universal cover R 2 , with slope 
a = U2/UJ1, assuming that uj\ 7^ 0. The main idea is to perform a change of 
basis enlarging the region of R 2 around the orbits, and a scaling of time. 

In the case a > 1, we choose the new basis {(1, 1), (0, 1)}, that decreases 
("D") the slope of the transformed vector: a' = a — 1. On the other hand, 
if < a < 1, one simply swaps ("S") the coordinates, i.e. we choose 
{(0, 1), (1, 0)} thus inverting the slope: a' = 1/a. The change of basis 
matrices are then D = and S = [?q], respectively. Note also that 

D a = [„ §], a G Z and 5 _1 = S. Finally, if a < 0, we recover one of the 
previous cases using V = [~q °] = V" 1 . 

Recall that the matrices D, S and V generate GL(2, Z). This group acts 
on]R = MU{±oo} asaM (c + da) / {a + ba) with [ a c b d ] 6 GL(2,Z). Each of 
the matrices D, S and V correspond to a transformation on K (in fact their 
inverses since those are the direct transformations to the new basis): D and 
S as above and V : a' = —a. 

The continued fraction expansion of a S R, is given by 

1 

a = a H , 

ai H ■ 

a 2 + ■ • • 

or simply a = [00,01,02? ••• ]j with ao G Z and a n £ Z + , n > 1. It is 
an infinite sequence if and only if a is irrational [|llj. We assume in the 
following that that is case and a > 0. In these conditions, the continued 
fraction expansion is a one-to-one map. 

Let [x] = max{A: £ Z : k < x} be the integer part of x S R and {x} = 
x — [x] the fractional one. The Gauss map is given by 

G: x i->- |i| , 2>0. 

Writing x n = G(x n -i), n > 1, with xo = {a}, one obtains the above 
coefficients a n from the recurrence relation ao = [a], a n = [l/x n -i], n > 1. 
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Note that a n = a n + x n , n > 0, where a n = [a n , a n+ i, . . . ], hence a n+ i = 
l/{a n } = l/x n . Define also 

n n 1 

^=n^=n— -■ (d 

1=0 2=0 

If Xk > 7" 1 for some fc = 0, ...,n — 1, where 7 = 1+ 2 V ^' , then, letting 
m = x^ 1 — Xfc+i > 1, x^Xk+i = 1 — mxk < 1 — < 1 — 7 _1 = 7 -2 . On the 
other hand, if Xk < 7 _1 , then x^+i > 7 _1 , and the product is again less or 
equal to 7~ 2 . Hence, 

f3 n <J- n , J^<J~ {n - j) , 0<j<n. (2) 
Pj-l 

The quadratic irrationals (roots of a quadratic polynomial over Q) have 
eventually periodic continued fraction expansion pi| , e.g. 7 = [1, 1, 1, ... ] 
and V2 = [1,2,2,...]. 

The action of D a ° on a is the elimination of the first coefficient ao of the 
continued fraction expansion, [00, a\, 02, • • • ] 1— > [0, ai, 02, . . . ]. The trans- 
formation S acts on sequences of the form [0, a±, 02, . . . ] giving [a%,a2, ■ ■ ■]■ 
The composition is simply the shift of the coefficients. 

We are looking for a specific type of linear coordinate changes of the 
torus. That is, it has to correspond to a matrix capable of improving the 
analyticity, as will be seen in the following sections (see also fl^l ). So, it 
has to be a hyperbolic matrix in GL(2,Z), with an unstable direction close 
enough to the subspace spanned by u. It turns out that we can choose the 
above shift of coefficients. The complete linear transformation is then 

T a = D a S= [ID e GL(2,Z), 

where a = [a]. If a > 1, this is a hyperbolic matrix with two real eigenvalues: 
one inside and the other outside the unit circle. More specifically, 

, a + vV + 4 

A = >7 

and — 1/A are the eigenvalues. The corresponding unstable and stable eigen- 
vectors are (1, A) and (1,— 1/A). 

Now, it is useful to recall some properties of the continued fraction ex- 
pansion, since we have associated it with the change of basis to be used in 
the construction of the renormalisation. We can write the above continued 
fraction expansion procedure in terms of the "convergent" matrices 



Pn 



7n-l Pn-l 
Qn Pn 



n > 0, and P_i 



1 
1 
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using the recurrence formula P n +i = T n+ iP n , where T n = T an , and 

Pn+1 = a n+ ip n +p n -l, 

Qn+i = a n +iq n + Qn-i, (3) 

fln-2 = dnPn-l + Pn, 

assuming = 1. The last equality is simply derived from rewriting P n -2 = 
OnA»-i- Notice that T* +1 = and P n = T n ---T , where T* 

stands for the transpose matrix of T. 

The continued fraction expansion can be obtained in only one step by 
a change of basis to {(q n -i,p n -i), (q n ,Pn)}, with the corresponding matrix 
being P*. Thus, (l,a n+ i) = (q n - X ~ aq n )P*(l,a). So, 

« = : , ( 4 ) 

q n -\ + q n a n +i 

which implies x n = -(q n a - p n )/(q n -\ - p n -i) and, by @), 

Pn = {-l) n {aq n -p n ). 
Again, from (^) using (n + 1) instead of n, \f3 n \ = l/(q n+ i + q n x n+ i). Hence, 

' < p n < — . (5) 



2^n+l qn+1 

It is easy to see that det P n = (— 1)" +1 and each line of the matrices 
P n G GL(2,Z) gives the convergents p n /qn = [a-o, ■ ■ ■ , a n] -> a, as n -> 
oo. This sequence verifies the best diophantine approximation property, 
\p n — aq n \ < \p — cm? |; for any integer p and < q < q n - Theorem 182). 
Also, if |q —p/q\ < l/2g 2 , then p/q is a convergent of a (|| - Theorem 
184). Using the fact that | det P n \ = 1, it is clear that the convergents verify 
\p n /q n -Pn+i/q n +\\ = l/q n q n +i- Then, 



Pn 

a 



< • (6) 

qnqn+l 



Denote the product of the first (n + 1) continued fraction coefficients 
of an irrational a by A n = nr=o a *- Rewriting @ one obtains q n +i/q n = 
o-n+i + V(<2'n/<Zn-i) ■ By induction it can be proved that the following is true: 

n 

A n <q n < A n JJ(1 + l/aiOi-x). 

8=1 

Similarly, we define by A n the product of the (n + 1) first numbers ce, = 
[di, aj + 2, • • • ]: 



in=n 
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This is related trivially to [5 n by the formula A n+ \ = ao/f3 n . So, from (||), 
K-i = a oA~ l + A~l v and, from (|), 

a q n < A n < 2a q n . (7) 

Also, using @, 

A n >a 0l n -\ ^>7"- J , 0<i<n. (8) 

Aj-! 

3 Space of Poincare Analytic Vector Fields 

Let (j)t be the flow generated by a continuous vector field X on T 2 , 6 = X(9), 
6 £ T 2 , and &t its lift to M 2 . For some norm || • || in R 2 , we say that 
wx(8o) = limt-^oo ^t(^o)/||^i(^o)[| is the winding ratio of X for the orbit of 
6o G T 2 , if the limit exists and lim^oo ||^t(^o)|| = 00 • Otherwise, if $t(0o) 
is bounded we put wx[8o) = 0. On the other hand, if the limit does not 
exist or if ||<3?j(#o)|| is unbounded but does not tend to infinity, we do not 
define the winding ratio. 

In the two-dimensional case, the set of winding ratios of vector fields X 
on T 2 is always a subset of {w, 0, — w}, for some normalised vector w (cf. 
e.g. |l|, |3|). Such a set is called winding set, and will be denoted by w(X). 

For any vector field X, the existence of an equilibrium is equivalent to 
S w(X), because it corresponds to a bounded orbit. If the slope of w S 
w(X) is an irrational number a, then w(X) = {w} (see || and references 
therein) . Fixed-point-free flows on T 2 with only one winding ratio are called 
Poincare flows. 

A closed cross section or transversal to a vector field on a surface is a 
simple closed curve, such that the vector field is nowhere tangent to the 
curve, and intersects all the orbits. A flow on T 2 has a transversal if and 
only if it is a Poincare flow. In this case, considering the return map on 
the transversal, it is possible to reduce the flow to a diffeomorphism of the 
circle. Therefore, there is an equivalence between these systems and the 
results are directly related. If w is rational (more precisely its slope), there 
is a periodic orbit of winding ratio w and every orbit is asymptotic to such. 
When w is irrational and the flow is C 2 , then all the orbits are dense on 
T 2 , and the flow is topologically equivalent to a constant flow with the same 
winding ratio || [l5|] . 

We briefly recall the space of vector fields used in jjTJ] that we will also 
study here for the two-dimensional case. 

Let r > and the lift of T 2 to a complex neighbourhood of M 2 : 

V(r) = {#eC 2 : ||Im0|| < ^} , 
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where || • || is the ^i-norm on C 2 . That is, \\z\\ = \z\\ + \z2\, with | • | the usual 
norm on C. We also denote the inner product as z ■ z' = ziz[ + Z2z' 2 . 

We will be dealing with analytic vector fields of the form X = u + /, 
where u G M 2 and with analytic functions /: T>(r) — ► C 2 , 27r-periodic in 
each variable 0j, represented in Fourier series as 

f(0) = £ / fc e™, A G C 2 . 
fcez 2 

More specifically, / belongs to one of the spaces (A(r), \\-\\ r ) and (A'(r), \\-\\' r ) 
whose elements are such that the respective norms 

11/11, = E H^ll er||/£|1 and ll/Hr = E ( X + 27r ll fc ID II^H erW 
fcez 2 fcez 2 

are finite. Both are Banach spaces. 

Some of the above vector fields X = uj + / generate Poincare flows on 
T 2 . Note that, if ||/||', < ||u;|| there is no equilibria. We will be interested in 
Poincare flows, especially with "irrational" winding ratio. 

4 One-step Renormalisation Operator 

In this section fix a vector uj = ((^1,(^2) £ R 2 \ {0} with u)\ 7^ and slope 
ol = W2/0J1 > 1, and denote a = [a]. As in JT0[| , we separate the terms of a 
vector field X into the following: 

Definition 4.1 For a > and i(j G C 2 , we define the far from resonance 
terms with respect to tp to be the ones whose indices are in 

I-{i>) = {ke1 2 :\i>-k\>o-\\k\\}. 

Similarly, the resonant terms are = Z 2 \ The respective pro- 

jections and on the space of vector fields are defined as 

[ltmX(9)= E X k e M . 

Let uj' = {a}~ 1 T~ 1 uj, with slope a! = {a} -1 that corresponds to the 
shift of the continued fraction expansion of a, so that ||u/||/||a;|| = (1 + 
a')/(l + a). Denote also by CJ = uj' /(uj' ■ uj') the normalised vector in the 
subspace spanned by uj' , i.e. Cj' ■ uj' = 1. We denote by K(X) the spatial 
average of a vector field, E(X) = L, 2 d(9, where d6 is the normalised 
Lebesgue measure. 

Proposition 4.2 There exists < k < 1 and a > suc/i i/iai, if0<p'<p 
with up < p' , any X G I+(u;).4,(//) has an analytic extension on T a T>(p). 
The linear map T a (X) = T" 1 X o T a from 1+ (ui)A(p') to A'(p) is compact 
with norm \\T a \\ < §ira/(p' — up) + 3a. 
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The proof of this proposition is contained in Section [|. As it will be seen, 
it is the existence of k < 1 such that I+(oj) C {k eZ 2 : < k||A;||} that 

guarantees the above analyticity improvement. 



The following theorem (proved in [12] using a "homotopy method") 
states the existence of a nonlinear change of coordinates U isotopic to the 
identity, eliminating all the far from resonance terms I~(u') = I~,,(u)' /a') 
of a vector field X in a neighbourhood of a non-zero constant vector field. 



Theorem 4.3 (gj) Let < p' < p, if) E C \ {0} and < a < \\if)\\. If X 
is in the open ball B C A'(p) centred at tf) with radius 

A y/6-2 f p- f/ 3-y/E a \ 

-a mm 



12 I 4?r ' 6 

there is a diffeomorphism U : V(p') — > isotopic to the identity, satisfy- 
ing 

[I~ (^)] (DU)~ X X o U = 0, and [7 = Id »/ [I" = 0. 

T/ie map U^: B ^ \^{il))]A{p') given by X i— > (Df/) -1 ^ o U is analytic, 
and the derivative at tp is Moreover, 



||W v ,(X)-^IU<2ri + max|^(3-v / 6),6(v / 6 + 2)M|^ ||x 



The definition below uses a different order of the linear and nonlinear 



maps, compared with the approach of [12]. There is no fundamental dif- 
ference, only a technical adjustment to allow us to evaluate the size of the 
iterates right after performing^. 

Definition 4.4 A map TZ^ acting on vector fields X is called a one-step 
renormalisation operator if it is of the form 



1 

Co' ■ E(X') 



K„{X)= - ^ 1', with X' =U UJ , /a ,oT a (X). 



The nonlinear transformation U u i / a i that eliminates the far from resonance 



terms I a (w'), is defined according to Theorem O with respect to u'/a' 



The linear map T a is given by Proposition 4.2. 



Since H, w {uj) = u>' , fixed points of IZ^ are possible for constant vectors 
to such that u' = uo. Those correspond to the ones with slope a having a 
constant continued fraction expansion, i.e. of the form [a, a,...], a > 1. 
Other periodic points of 1Z W can be obtained for a with a periodic con- 
tinued fraction expansion, which, for period p, is of the form [ao, . . . ,a p ], 
ao, . . . ,a p > 1. Eventually periodic points appear for slopes of vectors with 
eventually periodic continued fraction expansion (see |I2|| ). 
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Proposition 4.5 Given p' > and a > sufficiently small, we can find 
C > sitc/i £/iai the one-step renormalisation operator TZ^ is a well-defined 
analytic map from an open ball B of I^(cu)A(p') centred at u with radius 
Q = Ca 2 /{a\\J\\), toI+(uj')A(p'). In addition, ||7^(u;+/)-u/|| p / < \\f\\p>/( 
and \\n„(u + /)- J - DUM f\\ p , < [C(C - ll/IU)]~ lim|2 



Proof: Suppose p > p' and k < 1 such that up < p'. We will always make 
these choices such that p' — up is bounded away from zero. Then, by Proposi- 



tion H| T a (B) C B', with B' = {X G A'(p): ||X - w'/a'l^ < CCV/llw'll}, 
where C" is a positive constant. Then, C has to be chosen sufficiently 
small so that the transformation U^i i a i is valid in a neighbourhood of u'/a' 
containing B' , according to Theorem [D^. The same theorem states that 
K'/a'(B') C B" where B" = {X G 1+,, (u//a').40c/) : \\X - u'/a'\\ p , < 
CC'C"a/a'}, where C" > is an independent constant. 

Let F be the complex- valued continuous functional F{X') = u' ■ E(JT'), 
with the domain being a neighbourhood B'" C «4.(/f/) of uj'/a' such that 
£?" C B'" and J F(B ,// ) c{z£C: |a'z - 1| < 1/2}. Writing f = X'- u'/a', 
it is enough to have ||/'|| p , < £(4y/2-5)/2 < l/(2a'||cu / ||) with J = 1(1, cJ), 
which is satisfied for a sufficiently small C. Hence [&' • E(X')] -1 is bounded 
and analytic in B" and, for X' G B" , \\X'/(u' - EX') — u>'\\ p > < CC'C"C"'a, 



for some independent scalar C > 0. This, Theorem 4.3 and Proposition 



4.2 prove the first part of the claim. The first estimate also follows in the 
same way. 

To show the estimate on the second order remaining of the Taylor ex- 
pansion of 1Z W around u, consider g: z ^> 1Z w {ijJ + zf) — cj', f G B. This 
is an analytic map on an open ball in C with radius r = C/ II /lip' > 1- So, 
Cauchy's formula gives 

\\g(l)-g(0)-g'(0)\\ p , < ^ I ,^"\, dz 

2vr / [je[=r \z\z - 1)| 

< / 1 T\ SU P 11^(^ + -^'llp'- 

□ 



5 Diophantine Numbers and GL(2, Z) 

Let us review the definition of diophantine numbers and some of its proper- 
ties, and introduce the set of vectors of the plane with diophantine slope. 

Definition 5.1 For an irrational a and (3 > 0, we say that a is a diophantine 
number of order (3 if there is a constant C > such that 
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for any p/q S Q. We denote the set of all the diophantine numbers of order 
(3 by DC '{(3). All non-zero vectors uj in R 2 with slope a £ DC {(3) form the 
set denoted also by DC '(/?), which is GL(2, Z)-invariant. This is equivalent 
to say that we can find C > such that \u> ■ k\ > C\\k\\-( 1+ P\ k £ 1? \ {0}. 

If a 6 DC ((3) and using (0), there exists K > 1 such that 



1 

> 



Pn 

a 

q n 



K 



-i 



>~2+^ ri>0, 

q n 



where p n , q n are the convergents of a. This yields, together with (|3|), (|B|) and 
(0), that the diophantine condition can be written in the following ways: 

q n+ i < Kq n +I3 , 
a n+ i < Kq%, 



in+l 



< 2KAi +P , n > 0. 



Finally, for j3 > 0, one obtains lower and upper bounds for the conver- 
gents F n < q n < ^[(i+/3) n -i]//3, a nd for the coefficients a n+1 < K^+W, 
where F n is the Fibonacci sequence given by the recurrence formula F n+ \ = 
F n + Fn-i (i- e - with q n = F n and a n = 1). The sequence F n has an expo- 
nentially growing solution of the form 7™. All the constant type numbers, 
i.e. the ones inside DC(0), have an exponential bound on the sequence q n 
similar to F n , and a n < K. 

The set of the vector fields X close enough to E(X) with winding ratio 
w £ DC {(3) is Tv^-invariant for any uj G DC {(3). This is easily seen by 
recalling that if X' is the new vector field after a coordinate transformation 
isotopic to T € GL(2,Z), then the winding ratio of X' is Tu;/||Tu>||. In 
particular, DC '{13) is T a -invariant, and T a is the linear map involved in IZ^, 
modulo a constant rescaling. 

The constant type numbers DC(0) have zero Lebesgue measure, but 
the Roth numbers (see e.g. [||, ||): fl^oDC^), have measure one, which 
implies the full measure of any set DC '(/?), [3 > 0. Note that if Pi < P2 then 
DC{fr) C DC(p 2 ). 



6 The Limit Set of the Renormalisation 

Let wo £ K 2 \ {0} with non-zero slope «o = [o-o, a\, . . .]. For an initial vector 
field Xq, define the sequence of its images under renormalisation as X n+ \ = 
7Zoj„(X n ), whenever possible for n > 0. We denote also uj n+ i = 7£ Wn (u> n ). 
Notice that, if w(X n ) = w(tu n ), then w{X n+ i) = w{uj n+ i). 

The dynamics of uj n with slope equal to a n = [a n ,a n+ i, . . .], is given 
by the shift of coefficients, related to the Gauss map. If ao is a rational 
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number, then its continued fraction expansion is finite. Therefore, it is not 
possible to perform the iteration of this procedure for an infinite number of 
times. 

We want to include in the domain of this renormalisation iteration, vec- 
tors ujq with any slope and vector fields Xq with also far from resonance 
terms. So, initially, we allow TZ Wo to be an "adjustment" procedure, using 
the change of basis V or S and U. This unique transient step does not 
influence the qualitative result that follows. 

Theorem 6.1 Let p' > 0. For any vector ujq G U^<iDC(/3), there exists 
an open neighbourhood B of ujq in A(p') such that: if Xq £ B and w(Xq) = 
w(u>o), then 

\\X n -uj n \\ p , < K6 n , n>0, 
where K > and < 6 < 1. 

The proof of this theorem is included in Section |8[ For "irrational wind- 



ing ratios" not in the condition of Theorem 5.1, the domain of the renor- 
malisation shrinks faster than the rate of convergence of the iteration. Both 
of these behaviours are determined by the size of the coefficients a n . If their 
growth is controlled as for the diophantine case, the renormalisation reduces 
the perturbation enough to keep the iterates inside the domain. 

A corollary of the above theorem relates to the main result in [1^] for the 
two dimensional case. In fact, quadratic irrationals have an eventually peri- 
odic continued fraction expansion and are inside DC(0). It is then enough 
to choose a renormalisation operator 1Z associated to such uj, with the block 
of matrices T corresponding to the periodic string. In this way, u is an 
eigenvector of T and a fixed point of 1Z. Theorem |6.l| claims the existence 
of an invariant set contracting towards uj, that includes all vector fields in 
some neighbourhood with same winding ratio. To show that that is a sub- 
manifold, and there is a one-parameter family corresponding to the unstable 
direction, we refer to the details of the proof (namely the spectral properties 
of the derivative of 1Z at uj). The fixed point is therefore hyperbolic. 



Another application of Theorem 6.1 is a result for the Poincare map on 



the circle, in the same spirit as the one proved originally by Arnol'd |2| in the 
analytic case: any real analytic circle map with diophantine rotation num- 
ber a, inside some neighbourhood of the rotation R a : mi + a mod 1, 
is analytically conjugate to R a , for a norm induced by the one used here 



for vector fields (cf. Moser [16] for the differentiable case, and Herman ||] 
and Yoccoz [18] for more general and global results). We have not deter- 
mined the optimal condition obtained by Yoccoz |ll| in the context of the 
linearisation of germs of analytic diffeomorphisms and the local theorem on 
analytic conjugacy of circle diffeomorphisms. That condition states that it 
is necessary and sufficient that the set of winding ratios for which the claim 
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in Theorem 3.1 is expected to be valid strictly contains Up>oDC(f3), with 
slope verifying the Brjuno condition : Yl n >o ^ S(ln+l) / Qn < +°o. 

7 The Sets I+(u) and J£ 

Here it is enough to consider lu = (wi,^), ^1,^2 > 0, with slope a = 
U2/UJ1 > 1, and a = [a]. 

Definition 7.1 For k > define the set 

I K a ={kel?: \\T*k\\<K\\k\\}, 

and the projection I* on A(r), r > 0, given by 

r a X{6) = X k e 2mke , X e A{r). 

Let < a < uJi, i = 1,2, and 1/2 < k < 1. From the respective 
definitions one can rewrite the form of the cones I^{oj) and I£, closed subsets 
of Z 2 . In fact, is bounded by the lines k^ = mk\ and k% = lk±, and 

by fc2 = sfci and &2 = r&x, where 

_<£l^£ / = _<£!+£ and 1-K ±±K ( 1Q ) 

0J2+O"' OJ2— o" a— k V / 



Note that Z < m < and r<s<0ifK>(a + l 



1-1 



Lemma 7.2 If K> 1 — min{a(wi — a), 2(^2 — c) — a(wi + a)}, then 

Ia(oS) C /„. Moreover, considering uj = (l,a) u>i£Zi a > 1 and cr < 1/3, i< 
zs sufficient to have k > 1 — rj where n = (1 — 3<r)/3 > 0. 

Proof: It is enough to check the values of k for which we have r < I < m < 
s, as given in (10). The result follows from a simple calculation, considering 



the conditions on a imposed before. □ 



7.1 Analyticity Improvement in J* 

In order to prove that Proposition |4.2| holds, we need to show that X o T a 
is analytic in 2?(p) and has bounded derivative. From the result of Lemma 
|fj it is enough to notice that, for X(0) = £ fke 2nik ' e , 

IS 

\\XoT a \\ p < £ \\h\\e p ™ < £ llAlk^ 11 * 11 ^ 11 * 11 < 11*11,, 

and 

\\D(XoT a )\\ p <2^E ||T;fc|| e - <5 ll fc ll||/ (fc || e ^ +<5 -^ll fc ll e "'ll fe ll < y*pf[|,, 



11 



by choosing < 5 < p' — np, and making use of the relation sup t>0 te & < 
l/£for£ > 0. These bounds imply that \\T a {X)\\' p < (l+27TK/5)||T a ~ 1 || ||X|| p /. 
We choose e.g. 5 = n(p' — np). 

The above for T>(p) is true also for T>(r), r > p such that p' > rn. 
Therefore, T n = I o J, where J : I^A(p') — > A'(r) is bounded as T a above, 
and the inclusion map I: A'(r) — > A'(p), I(X) = X\ v ^ is compact. 



8 Convergence Result 

Here we include the proof of Theorem |6.1| , the convergence of the iterates 
X n towards the projected space KA(p') = C 2 of constant vector fields. 
More precisely, we show that the orbit of Xq approximates exponentially 
the orbit of lu . For all steps of the iteration, we choose the same k and p 



from Proposition 4.2. 



Considering the n-th step, the slope of u> n is given by a n > 1, and we 
denote its integer part by a n = [%]. While a n > 1 for n > 1, we can have 
ao > 0, but this can be regarded as a transient step (using S) which does 
not carry any problems into the following analysis. Also, if ao is negative, 
the change of basis V can be used initially, returning to the case of positive 
slope. Therefore, we will assume here a n > 1, n > 0. Again, the same 
can be said about starting within a space of vector fields with more than 
resonant terms. After an initial transformation, the iterates will always be 
inside a space restricted to vector fields with no far from resonance terms. 



Let a > as given in Proposition 4.2. To simplify the notation replace 



1Zu n simply by lZ n and use the projection P n of A(p') over uj n : 

where v = v/(v • v), for every v G M 2 \ {0}. According to Proposition 4.5, 



the domain of each 1Z n is the ball B n C A(p') around u n with radius Q n . 
The derivative at uj n is the linear operator given by: 

DKn(un) = (I-P n+1 E)£ n , 

with C n = a n+ il+ (uj n+ i) o T an . This is a compact operator since I + is 
bounded and T an is compact. 

The following proposition shows that there is a "super" exponentially 
shrinking of the non-constant terms, the image of (I — E), corresponding to 
the zero eigenvalue of the stable subspace. Recall that A n = YYi=o a i an d 
aoq n < A n < 2aoq n , with q n growing at least exponentially with n. Also, 
let I_i = 1. 

Proposition 8.1 Suppose w = i(l,a ) G DC((3), j3 > 0, t / 0. If a < 

|(1 + 47 3 ) -1 , then there exists c\,C2 > such that 

\\C n o ■ ■ ■ oCj (I - IE) || < cie- C2A ^||£ n o • • • o C j+1 (I - IE) || , 
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for j = 0, . . . , n with n > 0, and where 



A 2 ^ 



J 



Proof: Consider a vector uj with slope a, a = [a], the matrix T a and 
A = (a + Va 2 + 4)/2 > a. Let £1 be an orthogonal vector to u, e.g. = 
(1, — 1/a). So, 

T ( T 1 uj = ^-uj' and T a fi = — fi 7 , (11) 
a' a 

where O' = (1, — 1/a'). The projections of £ over ^ an( ^ ^ are 

= (u • k) uj and (I — W w ) k = (Q ■ k)£l, respectively. 
Introducing the notation 1+ = I^{uj n ) \ {0}, n > 0, and the correspond- 
ing projection I+, define a subset of if as 

= {k £ /+: T*---T*ke I+ V j <i<n}, 

for j = 0, . . . , n when n > 0. This set includes all the Fourier modes in 
Ij~ that will be mapped into I^+i (the only relevant ones since C n includes 
the projection by the sequence of matrices T* ■ ■ -T* . Therefore, for 

every k in Vj n , we can write \u> n +i ' T* . . .T* k\ < o~\\T* . . . T* k\\. Notice 
that \u n+1 -T*...T*k\ = \Tj ...T n uj n+l ■ k\ and, from ([□]), 



Tj . . . T n LOn+l - 

It follows that 

n+l 

\u> r k\ J] ai <a\\T*...T*k\\, keV+ n . (12) 
i=j+i 

A formula for the calculation of T* . . . T? k, with k £ Vy~„, can be obtained 
by induction: 




T*...T*k= (ujj ■ k) 



n+l 
i=j+l 



n 



+ n n+1 £ (O m+ i ■T^uj m )\ n <*i EI (-ai)' 1 

m=j \ i=j+l J i=m+l 

n 

+n n+1 (Cij-k) n(-ai) -1 , 

i=j 

(13) 
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making use of the orthogonality between uji and Qj for any i, i.e. uji ■ = 0. 
That yields T* Ui = (Q i+r T* u>i) + -If cDj) with u; i+1 -T* = 
Tj • a)j = There is K > 1 such that the norm of T* . . . T* k as 

given in (|i~3|) is bounded by 



\T*...T*k\\< \ Uj -k\ 



n+l 

ll^n+l|| 11 «i 



+2 £ ifw-r^j n si n <h l 

m=j \i=j+l I i=m+l 



+2K[U 



. -l 



a- 



We now need an upper estimate of the summation in the above inequality. If 
a is a quadratic irrational (to and are orthogonal eigenvectors of T) , then 
O m _)_i -Ti^ u m = for every m, which is in agreement with On the other 
hand, for other values of a, notice that | £l m +i -T^ Cj m \ < 1 ' ~ 1 
In addition, 



Tl I TfX \ Ti Tl ~\- 1 Tl 

e n °i n ar 1 < <*nli n «< e a^ 1 

m=j \ i=j+l I i=m+l i=j+l m=j 

n+l 

<7 3 "n+i 11 an, 
i=j+l 

because, from the second inequality in (||), we have that 



7 3 - 



The hypothesis and the formula oj ti = £(l,a n ), where ujq = £(l,ao), imply 
that T n a < 1/2 with T n = \\u n \\ +4>y 3 /(a n £). Then, using (g2|) it follows 
that 



2Ka\\k\\ 
\ujj ■ k\ < 

i - r n+ i<7 



n< 



aiai + i) 



(14) 



Given j > 0, as u>q 6 DC ((3), there is a constant Co > such that 



\ Uj ■ k\ = \u ■ tt 1 • • • t* = i k\f[ ai > f° n ;:i": (is) 

ll J ■ • • K W 



i=l 
72 I 



Consider the supremum norm on 

to || • || since \\\ ■ \\ < || • ||oo < || • II- Then, \\T^~ l . . . T*Z\ k\\oo = \\Pj~-i_ 
max{pj_i + Pj-2, Qj-i + Qj-2} < bAj—i ||fc||oo, for some b > 0, implies that 



! Halloo = rnaxj = i.2 \ki\. It is equivalent 

* — 1 rr*-l 7 II 

00 — II ■* i—-\ ro \\oo 



\T*- 1 ...T*zlk\\ < 2bA j ^ 1 \\k\\. 



(16) 
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Now, it follows from inequality ( |T4| ) together with (|l5|) and ( |l6| ) that, 
for all k € V£ n , 

with C = C /(2 3+ Pa b 1+ PK). 

Let the operator Vj n - A{r) — > «4(p')> r > p'j be a projection over the 
indices in Vj n together with an analytic inclusion. That is, 

where ||V^J| < exp [— C 1 ^ 2+ ^(r — p')Aj jTl \ . It is possible to include this 
transformation in the calculation of the sequence 

C n o • • • o Cj (I - E) = Cn o • • • o C j+ i (I - E) Yj~ n £j, 

< j < n, where £j : ^4(p') — > ^4(r) with p' > /cr, is £j followed by an 
analytic extension. The norm of Cj can be estimated above by ctj+ictj up 
to the product with a constant. Therefore, 

\\C n o ■ ■ ■ oCj (I - E)|| < ci exp(-c 2 A j , n )||£„ o • • • o C j+1 (I - E)|| 

with ci,C2 some positive constants independent either of j or n. □ 

To complete the description of the eigenspaces of DlZ n {oj n ), we look at 
the constant terms, the image of E. That is, we have to solve G n {f) = A/, 
with 

G n {f) = DKn{u> n )E (/) = (I — Pn+lK+lT^E (/). 

This operator can be re-written in the following way: 

-a n 1 
{a n }a n -{a n } 



Gn{J) 



fo 



l + {a n } 2 

As it has zero determinant its eigenvalues are and v n = tr(G n ), with 
\v n \ > 1 for any a n > 0. The respective eigenvectors are uj n and 
which determine the stable and unstable subspaces corresponding to the 
constant term. 

The constant part of each iterate X n = uj n +f n can be written as E(X n ) = 
w„ + P„oE (f n ) + (I - P n ) o E (/„). Notice that 

O E Un+l) = Pn+1 O G n {f n ) + 0„(||/X) 
= O ft ([|/n||J) J 

as we know from above that the image of G n is in (I — P n +i) C 2 . That is, 
its size can be estimated by the square of the norm of f n . The other term, 
(I — P n ) o E(/ n ), can be controlled in a different way, as it is possible to 
exclude some vector fields in B n that do not have winding ratio w n /||a; n ||. 



(17) 
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Lemma 8.2 Suppose that X E B n has winding ratio uj n /\\uj n \\. Then X 
belongs to the subset 

C n = {X£B n : ||(I-P n )oE(X)|| < \\(I-E)X\\ p ,}. 

Proof: A subset of vector fields D n C B n that do not cross the line spanned 
by ui n can be of the form: 

D n = {XeB n : \\X(9)-EX\\ < ||(I - P n ) o E (X)\\, Q E V{p')} . 

That is, the non-constant part of X E D n at each point 6 is less than the 
distance (given by the norm || • ||) between E (X) and its projection by P n 
over the subspace spanned by u n . The slopes of all the vectors X{9) are 
bigger than a n or always less than a n , never crossing that value (as for their 
respective winding ratios). Therefore, since we have that ||^(#) — EX|| < 
||(I — E)X|| p ' for every 9 E D(p'), the complementary set of D n in B n , 
contained in C n , includes all (but not only) vector fields with the same 
winding ratio as u; n . □ 

One can determine the non-constant part of each iterate X n = uj n + f n 
by the recurrence formula: 

(I - E) f n+1 = C n {l - E) f n + (I - E) O n (\\f n f pl ), n > 0. 

By induction, we obtain 

(I — E) f n+ i = C n o ■ ■ ■ o C (I - E) f 

+ £ C n o • • • o Cj (I - E) Oj-i (H/i-iHj) (18) 
j=i 

+(I - E) O n 

Note that the symbol Cn(||/n||p/) denotes the quadratic remaining of 1Z n 
as in Proposition which we write as F%. 

Lemma 8.3 If uq E DC((3), (3 < 1, there exists ct,C,t > such that, if 
H/olU < C 2 < I, we have \\f n \\ p , < CA~ T <( n ,n> 1. 

Proof: We prove this lemma by induction. From the formula f n = (I — 
E)/ n + P„, o E (/„) + (I - P n ) o E (f n ), and referring to (W7\) and Lemma K%, 



\\f n \\ p , <2||(I-E)/„|| p , +F- 



n—li 



where F 2 _ 1 = C„ii(Cn-i- ||/n-i||p') ^l/n-illjy according to Proposition |45. 
We write ( n = d / (a n a n+ i) for some constant d > 0. Note that C / A n < |Cn 
for a good choice of C and r independent of n, and by using the diophantine 
conditions. 
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Notice that e 1 < (v/t) v for any t > and v > 0. Therefore, for r > 



and < j < n, and the constants C\,C2 > given in Proposition 3.1: 



< Cl<7 7 



r(2 + /3) 



r(2+/3) 2^(1+0) 1/(14-/3) 



a: 



n+l 



(19) 



n+l 



where we have assumed a T < g(c2), with a function g(y) = [y/(r(2 + 
/ g))] r ( 2 +^)/(7ci), y > 0, for the last inequality to hold. 

If n = 1, provided that C < ±Co and 2C 2 (1 + 1/Co) < C/i{ < ±Ci, we 
have 

C 2 , C 1 



'Hip' 



< 2C 2 + 



< 1Z < TtCi- 



(o((o -C)~Al 2 

Assuming the result is true for n > 1, it remains to determine a bound of 
||/n+i||p'- So, we look at each of the terms of the formula (|l8|): 



E) /n+l||p' < o • • • o £ (I - E) /o|L' 



" I] ||Ai ° 

-F 2 



^(i-E)!!^ 2 



(20) 



The first two are estimated by the use of Proposition 8.1 through (fTsj) . 
Therefore, 



Co (I-E)/ |L< Cl Ce- 



< 



c 



7A 



n+l 



For the remaining terms we make use of the diophantine conditions in (^ 
with constants K, K and j3: a n+ \ < KAn and A n+ i < KA n + ^ . So, 

?2 2C 2 a\ a 2 +1 ^ 2C 2 K 2 a 2 if 



< 



< 



i 2 - " C 2 A% 



,j2 



T(2r-2-2/3)/(l+/3) 
^n+l 



< 



c 



^n+1 



for suitable C > and r > (2/3 + 2)/(l - /?). Now, 
Ej=i ll^n ° ■ 



„ „ rr(l+/3)+2/3+2 

^•(i-E)!!^ < 

71+1 7 — 1 



< 



+ 1 3- 

2C 2 K 2 \pn 7_t 7-t 
T2 - 2^7=1 



(21) 



7c'' 



n+l j— 1 5 



if r > (2/5 + 2 + t)/(l - /5) for some t > and /3 < 1. Since Ij^ can be 
bounded by times a constant, the sum in (^) can be bounded by 

some M > 0. Hence, (21) is estimated by C/{7A 1 n+1 ) as long as C is small. 
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Finally, if C, r and a are chosen accordingly to the various conditions 
described above, then 

2(C/7 + C/7 + C/7) + C/7 C I 

II /n+l II// < = < «Cn+l, 

as required. □ 
From the lower bound on A n in (||), ||/ n ||p' decreases with n at least 



geometrically like 7 Tn . That completes the proof of Theorem 6.1 
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